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1. The basic conceptions 

Let us consider the movement of particles (vehicles) on multilane 
road fragment. Let v be a velocity of regular movement deter- 
mined by a number of particles called slow vehicles that differs 
sufficiently from zero. Let d = d(v) be a dynamic distance, that 
includes one lane part of a road with the length that covers the 
length of vehicle and braking way (the size of discretization step 
by space coordinate) [1], fig. 1. 
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Figure 1: Discretization 

The dependence d(v ) can be approximated by the quadratic 
relation 

d(v) = c + c\v + c 2 v 

where Co is the length of vehicle base, c\ is a coefficient which is 
connected with regard of driver's reaction in case of unexpected 
traffic condition change, c^v 1 is an evaluation of braking way. 
For example let us consider Tanaka model, [2] 

d(v) = 5.7 + O.Uv + 0.0022*; 2 , (1) 
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d(v)(m), v (km/h). If in the equation (1) the velocity is mea- 
sured in m/sec, then Co = 5.7 m, c\ = 0.14 • 3.6 = 0.504 sec, 
c 2 = 0.0022 • (3.6) 2 = 0.0285 sec 2 /m. In this case we have 

d{v) = 5.7 + 0.504*; + 0.0285f 2 . (2) 

The coefficient c 2 depends on condition of the road covering. 
So, according [3], for the wet asphalt-concrete road covering the 
coefficient c 2 is in two times larger than for dry asphalt-concrete 
road covering, that is d(v) = 5. 7+0. 504^+0. 0b7v 2 , and for a road 
covered with ice d(v) = 5.7 + 0.504^ + 0.165i; 2 . Fig. 2 shows the 
dependence d(v) for different conditions of road covering. The 
graph of the function d(v) for different conditions of the road 
covering is shown on fig. 2. 
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Figure 2: Dynamic distance: di(y) - - dry asphalt; d 2 (v) — wet asphalt; 
ds(v) — road covered with ice 

Thus the immediate location of a vehicle on the road can be 
shown by the cell field, which for simplicity is synchronized in 
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Figure 3: The three-lane cell field 



Let T be magnitude of discrete time unit, that fixes flow 
state. If the "snap" does not change then the flow is supposed 
to be steady. The ratio of the number of occupied cells to their 
total number is called the regularity, r E [0; 1] . If r = 1 then 
we have the steady flow, an army column, which presents the 
movement of the column with the constant velocity v = const. 
It is a consequence of busy cells. If r < 1, then the individual 
transition of a vehicle to a front or neighbour cell diagonally is 
possible in a time unit T (fig. 4). 
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Figure 4: Manoeuvres on the road 



These transitions are caused by several reasons but one of the 
most important is that drivers want to drive their cars with a 
higher velocity. Let p = p(t) be a stochastic measure of individ- 
ual transitions at time unit t to the cell ahead when this cell is 
empty. Of course this value also depends on other characteristics 
of the flow and this is the subject of the further consideration. 
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Thus each car will do an attempt to move forward with prob- 
ability p{T) independent of the behavior of the other vehicles. If 
r ~ then there is no obstacle for such transitions as a rule, 
and if r is essentially greater than then the considered mea- 
sure depends on r. Let p(r, t) be this measure. It is clear that 
p(0, t) = p(T) and p(r, t) is an non-increasing function r G [0; 1], 
p(l,T) = 0. 

Let us evaluate the function p(r, t). Suppose that in the neigh- 
bourhood of the considered cell the states of the three neighbour 
cells that follow ahead are independent (fig. 5). 
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Figure 5: Manoeuvres from an inner cell 

Thus accurate to the coefficient p(T) the probability Pi(r; T) 
of a transition equals 

R 3 (r) = (1 - r) • 1 + r(2(l - r) 2 - (1 - r) 4 ) = 

= (1 - r)(l + 2r(l - r) - r(l - r) 3 ) = 

= (1 - r)(l + 2r - r 2 - r + 3r 2 - 3r 3 + r 4 ) = 

= (1 - r)(l + r + r 2 - 3r 3 + r 4 ) = 

— \ 4r 3 + 4r 4 r^. 

Similarly, for the case of inside lane on a multilane road (or on 
a two-lane road), we have 

R 2 (r) = 1 - 2r 2 + r 3 . 
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At last for the movement on a lane (or for canalized move- 
ment) we obtain 

Ri(r) = 1 — r. 

Thus 

pi(r, T)~p(T)(l-r). 
p 2 (r,T) ~p(T)(l-2r 2 + r 3 ); 
p 3 (r, T) ~ p(T)(l - 4r 3 + 4r 4 - r 5 ), 

(% 6). 
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Figure 6: Evaluation Ri(r) of dependence Pi(r;T)/p(T) 

In common case it can be considered that p(r; T) is a contin- 
uous function on r; (the addition of one busy cell cannot have 
any essential influence upon the mean velocity). 

At last in common case the considered function can be also 
depend on v. Really, every driver has his own knowledge of ve- 
locity. Therefore the additive component to expected velocity 
(vi > v) is compensated by individual attempts. For example 
according to the next scheme we have v + d(v)p(r, T,v)/T = v\. 
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2. DST-flows (Deterministic-stochastic traffic) 



The regular velocity can be also considered as the determinate 
component Vdet of the flow. Beside the collective motive of the 
particle behavior in the flow each particle also has its own in- 
tensions. Hence the flow can be presented as a composition of 
total (common, socialist) and private (liberal, individual) behav- 
ior. Let us consider a flow which consists of the particles with 
identical strategies of behavior (homogeneous dst-flow,). Then 
the flow velocity is the sum of the deterministic and stochastic 
components 

Vdst = Vdet + V st (3) 

where the stochastic transitions are independent and distributed 
equally. In this relation average meaning of individual transition 
is 

Vst=p(r,T,v)-^. 

Hence the average value and the dispersion of the velocity is 
accordingly equal to 

v = v + p(r,T,v)^, (4) 

v*=p(r,T,v)(l-p(r,T,v))(^) 2 . 
Let p be the density of traffic flow per lane. Then 

T 



Hence we have intensity per lane 



r d(v^) 
q~ pxv = -^(v + p(r, T, v)-=-) = 
a{v) 1 

rv rp(r,T,v) 

+ ™ • ( 6 ) 



d{v) T 
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The function (6) depends on three variables. As 
Pi( r , T, v) ~ p(0, T, v)Ri(r) 
we have (i is number of lanes) 

rv rp(0,T,v)Ri(r) 



d{v) T 
Suppose T — > 0. Reassume that 

p(o, r, u) 



T 

Then 



p(v). (7) 



= p(0, T, v)Ri(r)^- = p(v)Ri(r)d(v) 
and the intensity is equal to 

TV 

Qi(v,p(v),r) = -jj-j+rp(v)Ri(r). (8) 

According to (5) we receive 

Qi(p, v,p(v), p) = pv + pd(v)p(v)Ri(pd(v)). (9) 

Equation (9) generalizes the classical relation (main diagram) 
which is obtained for p(v) = 0. 
In this case 

Vi = v + p(v)Ri(r)d(v). 

3. The single- lane traffic "Regularity- Velocity- 
Rate" 

Let us consider the case of one lane, n = 1. We have v s t = 
p(v)(l — r)d(v) and the equation (6) for small T can be written 

as 

TV 

qi ^ -77-r + p(v)r(l -r) =vp + p(v)d(v)p(l - d(v)p), 
d{v) 
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where pd{y) = r < 1. For pd{y) > 1 the dynamic distance is not 
regarded and for this reason v diminishes. 
As d(v ) = - = Co + c\v + C2V 2 we have 

rv 

qi - Qi(r,v,p) = ■ ■ ^+pr(l-r), (10) 

C + C\V + c 2 v 2 

where coefficients Co, ci, C2 are assigned as in equation (1); p = 
p(v) 1/sec. 

The function q\ is defined on the rectangular < r < 1, 
< v < v max . Suppose p = 1. Let us represent the graph of the 
intensity function (fig. 7). 

For the same values of parameters the dependence of the 
individual velocity on r and v is shown on fig. 8. 
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Figure 7: Function q±(r,v, 1) vehicle/sec 



It can be seen from fig. 7 that the largest flow rate is reached 
when the values of basic variables r and v are near by (0,8; 15). 
In some neighborhood of point of maximum the trajectories of 
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Figure 8: Function vi(r,v, 1) m/sec 



intensity level lines are close of and are included to the consid- 
ered set of values. Thus the flow intensity can be still invariant 
when the parameters of the flow change. 

Comparing the behavior of graphs of the intensity and the 
velocity (fig. 7-8) we note that the conflict of the collective and 
individual purposes occurs. Just, if the velocity remains constant 
then the intensity changes, as the level lines intersect. 
In common case we have 
dqi v 



-^+p-2rp = 0, 
or Co + c\v + C2V 1 

dqi r rv(c\ + C2v) 



= 0, 



dv co + c\v + C2V 2 (co + c\v + C2V 2 ) 2 
i.e. cq + C\V + c — 2v 2 = v{c\ + 2c2v) and C2V 2 = Cq. Thus 



v = 



CO 
\C 2 
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4. Stability of the traffic and scattering of the funda- 
mental diagram 




Figure 9: Experimental data in field "density-intensity", [4] 



10 





Figure 10: Experimental data in field "density- velocity" 
with different intervals of averaging, 
the upper fragment is taken from [4], 
the below fragment is taken from [6] 
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The considered model allows to explain the appearance on 
the fundamental diagram of scattering domains of density where 
different values of intensity can correspond to the same value of 
density. The existance of scattered domains has been discovered 
on the basis of the experimental results, described in [4,5], fig. 
9-11. 

If the flow density p < /9 max /2 then the movement is stable 
and the traffic can be described well by a hydrodynamic model. 
In this case the regular movement r = 1 is a single stable state. 
The effect of random disturbances are short-lived and the traffic 
return fast to the previous regime. 

For sufficiently large densities, that are out of the range of 
the stability of the traffic, the appearances of disturbances result 
in that regularity of traffic has values r < 1 because of reaction 
of drivers on changes of the situation that gives a change of the 
dynamic dimension. The density becomes a random value. Dis- 
persion of the density becomes not equal to zero. The changes 
of rate r/d(v) result in that values of densities, which have been 
found in the experiments, has dispersion near the average value 
when the intensity of the traffic does not change. Thus the same 
value of the intensity, obtained in the experiments, can corre- 
spond to differents values of densities and therefore scattered 
domains on the fundamental diagram appear. 
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Figure 11: Fundamental diagram with scattering domains 

5. Multilane motion "Regularity — velocity — intensity" 

Similarly, for the case of two lanes 

q 2 = q 2 {r, v,p) = ■ ■ ^ + pr(l - 2r 2 + r 3 ). (11) 

we have the graphic dependences for the intensity and the ve- 
locity (fig. 12, 13). 
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Figure 12: Function q<i{r, v, 1) (vehicle/sec) 
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Figure 13. Function ^(r, u, 1) (m/sec) 
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Finally, for three lanes 



Q3 = qs(r,v,p) = 



rv 



+ pr{l - 4r 6 + 4r 4 - r 5 ) (12) 



Co + C\V + c 2 v 2 
we obtain the dependences shown on fig. 14, 15 
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Figure 14. q^(r, v, 1) (vehicle/sec) 
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Figure 15. ^(r, u, 1) (m/sec) 



15 



6. Collective and individual for the flows with a con- 
stant density 



For fixed density po = r/d(v) from (12) we obtain 



qi = qi(p,v,p) 



= vp + pd(v)p(l — d(v)p) = 

Vp + p(c + C\V + C 2 V 2 )p(l - p(c + C\V + c 2 v 2 )) 
= d~ 1 (r/v)p + pr(l — r). 



(13) 



The unequality r < 1 is equivalent to unequality pd(v ) < 1, 

i.e. 



Let us represent the evaluations of dependence (13) of inten- 
sity on velocity for the fixed density. The question is how the 
rate changes if the flow velocity changes abruptly during a small 
period of time and the density is constant. If the case of one 
lane is considered then for po = 0.01 vehicle/m we have the de- 
pendence shown in fig. 16. As for dependences, shown on fig. 
17-19, this dependence corresponds to the case of dry asphalt 
(eq. (2)): Co = 5.7 m/sec; c\ = 0.504 sec; c 2 = 0.0285 sec/m. 
Range of velocities below 25 m/sec, i.e 90 km/h, is considered. 



C + CiV + c 2 v 2 < - 

p 




(14) 
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Figure 16. qi(v,p) (vehicle/sec), p = 0.01 (vehicle/m) 
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Figure 17. Dependence q~i(v,p) (vehicle/sec), p = 0.05 (vehicle/m) 
The dependence for po = 0.05 is shown on fig. 14. On fig. 17 
there is the evident unstability in the neighbourhood of the right 
lower angle. For the case of the three-lane road we have depen- 
dences shown on fig. 18-19. 
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Figure 18. qz(v,p) (vehicle/sec), p = 0.01 (vehicle/m) 
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Figure 19. qz(v,p) (vehicle/sec), p = 0.05 (vehicle/m) 
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7. Influence of "blue lights" upon the flow rate 

The particles of the two types characterized by of the different 
functions of dynamic distances d\{y) and d,2(v) are considered. 
Suppose d,2(v) < d\{y) for each allowed v. If v is the regular 
velocity of the unmixed flow with zero individual velocity of the 
unmixed flow, then by mixing between two large particles small 
particle would emerge (percolation). Hence the flow velocity 
can be evaluated as v\{v) = d^ 1 (di(v) / (2)) . Let us suppose still 
that the number of such particles is rather small. That the flow 
intensity in the new conditions is 

where the function p(r, T) is defined as the function p(r, T) intro- 
duced in section 1; dynamic distance d(v) is calculated according 
equation (2). 

The value Vi(v), which satisfies equation d(v\(v)) = ^ that 
is equation 

2 c + civ + c 2 v 2 
CO + CiVi{V) + c 2 v 1 (v) = , 

exists only in the case 

civ + c 2 v 2 > cq. (15) 



Suppose also 



i.e. 



l-pd(v)/2 > 0, 
2 2 

CO + C\V + c 2 v < -, 



Cl 

v < h 



Cl ^ 2 _ Co 



2C2 \ 2C2 C2 C2p 
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Thus the equalities 



Cl , ( Ci 2 C 
^— + \ (75— ) + — < U < 

2c 2 \ 2c 2 c 2 



ci 



+ 



are true. 
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Figure 20. Dependence Q(p,v) 



Suppose that 



p(v)d(v)/2 = v — v\{y). 



Hence 



Qi(p, v ) = P ( v i( v ) + (v- v 1 (v))R l (pd(v)/2)) . 
For example, if i = 1 we have 

qi(p,v) = piyxiy) + {y - Vi(v))(l - pd 1 {v)/2)). 



As 



ci 



Vi(v) = ~- h 



, Ci C Ci v\ 



2c 2 V2c 2 ; V 2c 2 2c 2 " 2 
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we obtain the ratio of intensity of disturbed and undisturbed 
conditions 

= «iM + («-«iM)(i-MM/2) 

Relative variation of flow rate of slow particles in case of 
small number of fast special vehicles is shown on fig. 20. We 
still suppose that Co = 5.7 m/sec; c\ = 0.504 sec; C2 = 0.0285 
sec/m. 
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